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Abstract 

00 ■ 

■ From velocity-gated small-angle correlation functions the emission 

f"^ | chronology can be deduced for non-identical particles, if the emission is 

^\ . independent. This is not the case for non-identical particles that originate 

from two-body decay of fragments. Experimental results may contain 
, contributions from both independent emission and two-body decay, so 

care is needed in interpreting the velocity-gated correlation functions. It is 
shown that in some special cases, it is still possible to deduce the emission 
QJ , chronology, even if there is a contribution from two-body decay. 

o ■ 

2 ! 1 Introduction 

a ■ 



Understanding the different origins of emission and the emission time sequence 
of different particles is a major challenge of intermediate energy heavy ion col- 
lisions. For more than 20 years, two-particle intensity intcrfcromctry has been 
used as a probe, yielding a convoluted space-time information that, however, is 
difficult to disentangle El EI E] • This is because of the presence of multiple 
sources, and the competition of non-equilibrium emission processes with equi- 
librium relaxation modes, leading to a broad range of origins of the measured 
particles El El 11 M HH1 EH E21 El It is the purpose of this paper to discuss 
some of the assumptions made when inferring emission chronology from velocity 
gated experimental correlation functions. 

Particles emitted in nuclear collisions may interact with each other after the 
emission. This final-state interaction causes the population in phase-space to be 
altered, an effect that can be observed by the small angle two-particle correlation 
function. The proximity in space and time of the interacting particles influences 
the strength of the final-state interactions, and hence the size of the correlations 
seen in the correlation function C(q), where q = /Lt|pi/mi — p^/tn^l is the 
relative momentum and [i = toiTOi/(toi +1782) is the reduced mass. Correlation 
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functions constructed from experimental data therefore contains information 
of the space-time characteristics of the emitting source. Experimentally, the 
correlation function is constructed by dividing the coincidence yield, by the 
yield of non-correlated events, normalized to unity at large values of relative 
momentum, where no correlations are expected. Often, much of the space- 
time information contained in this 6-dimensional observable is lost, because of 
implicit experimental integrations over some of the dimensions in the relative 
and total momenta. When statistics is high enough, some of the information can 
be recovered by applying directional cuts ^3] an d total-momentum or energy 
gates |2 ■ Furthermore, for non-identical particle correlations, it is possible to 
apply particle-velocity gates and infer the emission chronology of the different 
particles QU EH El [13 d HJ . 

Velocity gated correlation functions of unlike particles is a very powerful 
method, as it is model independent. However, there are some critical assump- 
tions made, that may not always be fulfilled depending on the reaction scenario 
and selection of emitted particles. For the method to work, the particle veloci- 
ties should be obtained in the frame of the emitting source. The source velocity 
may not always be a well defined quantity, because of implicit integration over a 
limited range of impact parameters. In section|3|we investigate into some detail, 
consequences of uncertainties in the source velocity. Another crucial assumption 
for extracting emission chronology from the velocity gated correlation functions, 
is that the non-identical particles must be emitted independently by the excited 
source. If they instead originate from the two-body decay of a fragment, their 
respective velocity is determined solely by energy and momentum conservation, 
and they do not carry any useful information on the emission time. In section^] 
we investigate under which conditions the emission chronology may be inferred, 
when particles originate from both independent emission from a (large) source, 
and from two-body decay of fragments. 

2 Space-time characterization 

Normally, particles originating from a specific source are emitted at different 
rates during the time interval of emission, leading to a specific time distribution 
for the source. If a two-particle correlation function could be constructed from 
particles emitted from a single source, and if the spatial distribution would be 
known, then the shape of the correlation function would yield information on 
the shape of the time distribution 1 . Normally, the spatial distributions are not 
known. To interpret the experimental results, source models, that contain some 
assumption on the spatial and temporal distributions, are often used. The shape 
of these distributions can, to some extent, be varied by varying parameters of 
the model The best-fit parameters then represent the average emission point 
and average emission time, though it should be remembered that such average 

1 In Refs. 11 21 22 23 24) the relation between the shape of the correlation function and 
the shape of the source function is discussed. If the spatial part of the source is known, the 
shape relation in Refs. pill 12 II 22 23 24 is directly applicable to the time distribution. 
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values are model dependent. In the experimental data the situation is much 
more complex, since it is never possible to completely isolate one source from 
other sources present during the reaction. The contribution from several sources 
leads to a complex total time distribution. 

The term emission chronology is usually used to denote a difference in the 
average emission time between two particle types. The difference in the average 
emission times may be small compared with the width of the emission time 
distributions. A difference in the average emission times, extracted from exper- 
imental data, may have different origins, depending on where the particles, that 
are included in the observables, are coming from. Possible origins may be 

• A shift of two similar time distributions. This shift could be due to the 
nuclear interaction. For systems with an exotic isospin composition, the 
symmetry interaction could cause the emission time of neutrons and pro- 
tons to be different. 

• A different width of two, otherwise quite similar, time distributions. For 
an equilibrated source, differences in the Coulomb barrier for different 
particles, could lead to different widths of the time distributions. 

• Different relative weights of several sources. The abundance of emitted 
neutrons and protons from an intermediate velocity source, and non- 
equilibrated residues, could be different due to the isospin composition 
of the systems. 

When interpreting experimental data, it is important to be aware of different 
possible origins of different average emission times, in order to make the cor- 
rect conclusions. Different origins can be enhanced or suppressed by applying 
different gates and conditions on the observables. 

The emission chronology between two particle types (e.g. neutrons and pro- 
tons) can, under certain conditions, be determined from like-particle correlation 
functions (e.g. neutron-neutron and proton-proton). If it is valid to assume that 
both particle types are emitted from the same spatial region, a model fit to the 
experimental data, will yield an average emission time for each particle type. By 
comparing these average emission times, an emission chronology can be inferred 
|25| . The drawback of this method is that the results are sensitive to the as- 
sumption of emission from the same spatial region. Furthermore, the extracted 
average emission times are model dependent, since the average emission times 
depend on the shape of the (spatial and) temporal distributions assumed by the 
specific source model. 

3 Emission chronology from non- identical par- 
ticle correlations 

Model independent information on the emission chronology of two particle types 
(e.g. neutrons and protons) can be obtained from non-identical-particle correla- 
tion functions (e.g. neutron-proton). A technique was first suggested for charged 
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particle pairs, based on comparison of the velocity difference spectra with tra- 
jectory calculations [221123 EH- The technique was extended to any kind of 
interacting, non-identical particles, by applying energy or velocity gates, and 
proposed for particle pairs such as pd and np pw |2H], and K + K~ |3(J| . 

3.1 Particle-velocity-gated correlation functions 

The basic idea of velocity-gated correlation functions of non-identical particles 
[1611171 [T8| . is that, if there is an average time difference in the emission times of 
two particles types, there will also be a difference in the average distance between 
the particles, for particle pairs selected with the condition v± > v% as compared 
to the pairs selected with the complementary condition Vi < v 2 . It is obvious 
that in the class V\ > v 2 the average velocity of particle 1 will be higher than in 
the complementary class v\ < v 2 . This means that if particle 1 is emitted first, 
it will, with the condition v\ > v 2 , on average travel a larger distance before 
the second particle is emitted, than with the complementary condition. In this 
case the condition v\ > v 2 leads to on average larger distances (and weaker 
interactions) than the condition v\ < v 2 . 

The effect can be easily seen if one compares the correlation function C±, 
gated on pairs V\ > v 2l with the correlation function C 2 , gated on pairs v\ < v 2 . 
Assuming that particle 1 is on average emitted first, the ratio C\jC 2 will show a 
dip in the region of relative momentum where there is a correlation (attractive 
interaction) and a peak where there is an anti-correlation (repulsive interaction) . 
Furthermore, the ratio C\/C 2 will approach unity for both q — > (since the 
velocity difference of the two emitted particles is negligible) and q — > oo (since 
modifications of the two-particle phase space density arising from final state 
interactions are negligible). A single normalization constant, calculated from 
the non-gated correlation function, is utilized for both C\ and C 2 16 . Note 
that while the height of the non-gated and gated correlation functions depends 
on the normalization, and therefore is sensitive to the statistics, the ratio of 
the velocity-gated correlation is not sensitive to the normalization, since the 
common normalization constant cancels out. 

The exact location of the peak and/or dip in the ratio depends on the source 
and in particular on the origin of the difference in the average emission times. 
It should also be mentioned that it cannot be ruled out that the differences ob- 
served in velocity-gated correlation functions could have a spatial origin. How- 
ever, such a correlation would then mean that there is a correlation between the 
spatial region where the particles are emitted and the particle velocities, and 
that this correlation is different for the two particle types used in the correlation 
functions. For heavy-ion collisions at intermediate energies, such an explanation 
is clearly more unlike than a difference in the average emission times. 

3.2 Influence of the source velocity 

The main uncertainty in the method of velocity-gated correlation functions 
comes from the uncertainty in the source velocity. This uncertainty originates 
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mostly from an implicit impact parameter averaging in the experimental data. 
The range of selected impact parameters leads to a distribution of source veloc- 
ities. In addition, the measured particles in the pair could have been emitted 
from different sources with a relatively large difference in source velocity. By 
applying suitable conditions and gates on the experimental data, the fraction 
of such pairs should be low with respect to particles coming from the desired 
source. 

If the velocity of the assumed source is different from the real source velocity, 
the calculated particle velocities in the assumed source frame will contain some 
error, and it may happen that the magnitude of the two particle velocities is 
interchanged as compared to the real source frame. In this section we make an 
estimate of the fraction of pairs that end up in the wrong gate, depending on 
the differences between the real and assumed source velocities. 

(S) (S) 

Assume that the particles 1 and 2 are emitted with velocities and v 2 
from a source S. If the source has the velocity v« in the laboratory system, 
the particle velocities in the laboratory system become 

vf> = vf>+vW, i=l,2. (1) 

If the velocity ui L ' is used to calculate the particle velocities in the source 
system, we get 

= vf)-uW 

= v W + [vW-„W], , = 1,2. (2) 
The condition v[ S ^ > can a l so be expressed < [fi^] 2 — [v^] 2 , an d ^ 



is straightforward to show that 



+ 2[vW-uW].5. (3) 

The error of using the velocities Uj instead of the real velocities Vi can thus 
be estimated from Eq. ©. Noting that [v[ S) ] 2 - [v ( 2 S) ] 2 = [v[ S) + v { 2 S) ] ■ q//z, 
the expression shows that if a reasonably good source selection has been done, 
so that |vi L ^ — ui L ^| is small relative to |v^ + Vj |, the second term in the 
right hand side of Eq. J2J should be negligible. This is especially the case when 
the relative momentum increases, since |v 1 S ' ) + | increases with q for non- 
identical particles. 

As an example we have made a simple simulation for the case of one proton 
and one a-particle emitted from a source. We have assumed a Boltzmann distri- 
bution of the particle energies, with temperatures T(p) — 6 MeV and T(a) = 8 
MeV. 2 The particles have been assumed to be emitted isotropically. The dif- 
ference between the real source velocity and the assumed one has been taken 

2 We have assumed that the a-particle is emitted somewhat earlier than the proton 1191 . 
thus feeling a higher average temperature. However, the assumption of different temperatures 
does not influence the conclusions of this section. 
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Figure 1: Results from the simulation described in the text of subsection 13.21 
Panel (a) shows the fraction of pairs that are attributed to the wrong velocity 
gate due to the uncertainty in the assumed source velocity. Panel (b) shows the 
used distribution of source velocities. 



normally distributed, with random direction. The result is presented in Fig.^i, 
where it is seen that the number of pairs attributed to the wrong velocity gate, 
is around 5% for a realistic estimate of the uncertainty in the source velocity 
(see Fig. ^>). The error becomes larger if the source velocities have a larger 
spread, but even for an unrealistically large uncertainty in the source velocity, 
covering all possible velocities in the reaction, the fraction of wrongly attributed 
pairs is less than 25%. 

The example indicates that, for most cases, the error will be small, of the 
order of a few percent. Nonetheless, this should be explicitly verified in each 
application of the method of velocity-gated correlation functions. 



4 Two-body Decay from Fragments 
4.1 Background 

For those correlation functions characterized by final state interactions leading 
to resonances, the resonance peaks may a priori have two different origins (Hll 
E2EH1: 

1. From interactions between independently emitted particles from a (large) 
source 3 . 

3 We have here in mind that the two fragments interact and thereby change their momenta 
in such a way that the relative momentum region under the resonance peak is populated, but 
the particle do not form a fragment that then in turn undergo two-body decay. The latter 
process would not be distinguishable from two-body decay of fragments formed by other 
mechanisms in the reaction (corresponding to process gjin the list). 
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2. From processes where an unstable fragment formed in the reaction decays 
into the two measured particles (like in e.g. 8 Be -» a + a or 6 Li — > a + 
d). 

For the case when a fragment is decaying into two non-identical particles, the 
velocity of the two particles in the fragment rest frame is always such that the 
lighter particle gets a larger velocity than the heavier particle. This follows from 
energy and momentum conservation in the decay. 

If we could construct a velocity-gated correlation function with all particles 
coming from two-body fragment decays and with their velocities calculated in 
the fragment frame, we would find all the pairs in the gate where the lighter 
fragment has a higher velocity, and none in the complementary gate. In sub- 
section ED we show that even if the particle velocities are calculated in another 
frame, more than 50% of the events from fragment decay will still be attributed 
to the gate where the lighter fragment has the higher velocity. 

Any deviation from this behavior in experimental data would be due to 
particles from other sources than two-body fragment decays. Therefore, when 
in some cases it is observed that the gate where the heaviest particle has the 
largest velocity leads to a stronger correlation or anti-correlation than the com- 
plementary gate, it can reliably be concluded that this behavior is dominated 
by a mechanism different than two-body decay. In such cases the effect may 
be attributed to the interaction of independently emitted particles, and the 
velocity-gated correlation function can be used to obtain information on the 
time sequence of the independently emitted particles (see e.g. Ref. |19|). 



4.2 Kinematics 

Consider a fragment Ap with mass mp decaying into two particles A\ and A 2l 
with masses mi and m 2 , Ap — > A\ + A2 . Momentum conservation in fragment 
rest frame, 

pi = p 2 = q , (4) 

and energy conservation, 



m F c 2 + E* = y/( mi c 2 ) 2 + (pic) 2 + \J{m 2 c 2 ) 2 + (p 2 c) 2 

2 2 

^ 2 1 Pi 1 2 1 P2 /r\ 

* miC + 2^ +m2C + W (5) 



lead to 



(rriFC 2 + E*) 2 (m 2 c 4 — m^c 4 ) 2 m 2 c 4 + m 2 



c 



4 A(mpc 2 + E*) 2 2 

2m 1 c 2 m 2 c 2 (mp c 2 + E*) 2 2 



m\c 2 + 7712C 2 



2m\c m 2 c , (6) 



where q is the relative momentum of the two particles, and E* is the excitation 
energy of the decaying fragment. In the rest frame of the decaying fragment, 
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the velocities of the two emitted particles become 



i p ) 



= -q/mi 
= q/™2 , 



(7) 
(8) 



and the velocity of the lighter particle will always be larger than the velocity of 
the heavier particle. 



If the velocity of the fragment in the laboratory system is v F 
particle velocities in the laboratory system become 



' F 



I F ) - „w 



q/mi 



IF) 



v F + q/m 2 



then the 

(9) 
(10) 



In applications of velocity-gated correlation functions, the particle velocities 
are calculated in some source system (intermediate velocity, target or projectile 

residue source) with assumed velocity vi L ' . The particle velocities in the source 
system become 



I s ) 



i l ) 



i l ) 



— q/mi — v 



(L) 



(L) 

(L) (L) . i (L) 

v y s ' = v F + q/ m 2 - v ^ 



(11) 
(12) 



get 



The condition of the velocity gates, v\ > V2 (vi < v 2 ), can also be written 
> {v\ — v\ < 0) . Using the particle velocities in the source system, we 



/ (S)\2 / (S)n2 
(y\ ') - (v 2 ') 



m 2 — mi 
•m\m 2 



2e„ 



IS) 



(13) 



where e q is a unit vector directed along q, and = w y p ! — v\^' is the fragment 
velocity in the source frame. In a given decay, q, m\ and m 2 are determined by 

energy and momentum conservation, while the direction of q can be considered 

(s) 

to be isotropic. If we denote the angle between q and v^, by /3, the condition 
in Eq. 1)13(1 can be written 



I L ) ,I L ) 



(A 



(Sh2 



m 2 - mi (s) 

q - 2 v F 

m\m 2 



cos(P) 



(14) 



4.3 Results 

From Eq. i|14[l it is clear that in the source system, the light particle does not 
always have a larger velocity than the heavier. The condition depends on the 
source velocity, and on the angle between the relative momentum and the source 
velocity. 

In Fig. |21 we present the regions where the lighter or heavier particle has the 
largest velocity. When the fragment velocity (in the source system) is small, 
the velocity condition in the fragment frame still holds (i.e. the lighter particle 
has the higher velocity). But, as the difference between the fragment velocity 
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Figure 2: A graphical representation of the velocity condition in Eq. I|14(l for 
the decay 5 Li — > p + a. The gray area represents the region where — v 2 a < 0. 

and the source velocity becomes larger, a larger fraction of events will have 
the velocity condition interchanged. However, it is important to note that the 
number of events for which the original velocity condition holds is always more 
than 50% assuming that q is isotropically distributed. This means that if all 
particle pairs would come from two-body fragment decays, then the correlations 
in the gate where the lighter particle has the higher velocity (in whatever frame) 
will always be stronger than in the complementary gate. Thus if the contrary 
behavior is observed, the dominating origin of particles must be other than 
two-body decay. 

From En. ll4l it is clear that the regions in the source system (as in Fig. [2J 
where lighter or heavier particle has the largest velocity, depends on the relative 
momentum, q, in the decay. For two-body decay that yields large values of q, 
the region where the lighter particle has the largest velocity will be larger than 
for decays yielding smaller values of q. This dependence is the main reason for 
different results of the subsections I4.3.1I and I4.3.2I below. 

In experimental data, the situation is somewhat more complicated, because 
of limited angular coverage. The angular coverage may select certain regions 
in the cos(/3) - v F plane, that then could change the above conclusions. It 
is not possible to find a simple expression for the velocity condition including 
the angular coverage of a given experiment. Instead this has to be investigated 
numerically for each setup whenever velocity gates are used. In the next two 
sections we presents the results for two such investigations, and show how they 
could be implemented. 
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Figure 3: The experimental velocity-gated p a correlation function of Ref. 
The arrows in the figure indicates the positions of different fragment decays (see 
Ref. QU). 



4.3.1 Application to emission from projectile residue 

In Ref. velocity-gated small angle two-particle correlation functions were 
used for protons, deuterons, tritons and a-particles, from the E/A — 44 and 77 
MeV 40 Ar + 27 Al collisions, at very forward angles, with the aim of studying 
emission from the projectile-like residue (PLS). Fig.[3]shows velocity-gated p a- 
correlation functions for the E/A = 44 MeV reaction. One can note that in 
the region q ~ 54 MeV/c, corresponding to the two-body decay of 5 Li, the gate 
v a > v p shows larger correlations than the complementary gate v a < v p . One 
possible deduction from this observation (and the discussion from the previous 
section) would be that these particle must have a different origin than decay 
from 5 Li. However, before such a conclusion can be drawn, it must be ruled out 
that the very specific angular coverage of 0.7° < 6 < 7° in Ref. does not 
impose such constraints in the cos(/3) - plane of Fig. |2 that the velocity 
condition is interchanged. 

To investigate this issue, we have performed simple numerical simulations, 
where we have assumed that the proton and the a-particle come from decay 
of 5 Li, and that the 5 Li has been emitted from different sources. We require 
that the proton and the a-particle should come within the angular range of 
the detectors (0.7° < 8 < 7°), as well as be above the energy threshold of 35 A 
MeV, imposed in Ref. ^Hj- The results are summarized in Tabled The first row 
shows the results when we assume that the 5 Li has been emitted from a target- 
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Fragment 
source 


All 
pairs 

(PLS) ^ (PLS) 


0.7 < 

Detected 


: Blab < 7.0° 
4 > 35 MeV 

(PLS) ^ (PLS) 


TLS 
IS 
PLS 


54.3% 
58.6% 
76.9% 


0.00% 
0.03% 
2.27% 


84.7% 
82.8% 



Table 1: The fraction of pa pairs emitted in the two-body decay of 5 Li which 
have v p PLS ^ > v^ LS \ The last column gives the fraction of "detected" pairs. 



like residue source (TLS) of temperature 7 MeV. We deduce the kinetic energy 
of the 5 Li from a Boltzmann distribution, assume that the emission is isotropic 
in the TLS frame, and boost the 5 Li with the TLS velocity in the laboratory 
frame. The TLS velocity is taken normally distributed around the mean value 
of 0.02 c with a standard deviation of 0.01 c. The decay of 5 Li is assumed to be 
isotropic, while the magnitude of the proton and a velocities in the J Li frame is 
determined from the energy and momentum conservation in the decay. Finally 
the proton and a velocities are calculated in the frame of the projectile-like 
residue source (PLS), which is assumed to have the velocity 0.27 c along the 
beam direction. We have used 200 000 events in the simulations. The second 
column in Tabled shows in how many of all events, irrespectively of detection, 
that the proton has a larger velocity in the PLS-frame, than the a-particle. The 
third column shows in how many of the events that both particles fall within 
the angular range of the detectors, and above the imposed energy threshold. 
The fourth column shows in how many of the events in column three, that the 
proton has a higher velocity in the PLS-frame, than the a-particle. For the 
intermediate velocity source (IS) and the PLS we have assumed temperatures 
of 20 and 7 MeV, respectively, and source velocities normally distributed around 
the mean values 0.15 c and 0.27 c, respectively, with standard deviations of 0.02 
and 0.01. 

As seen in Table ^ the experimental filter in this case strengthens the ve- 
locity condition in the 5 Li frame, as compared to the unfiltered result (column 
two). This means that if the protons and a-particles in the region q ~ 54 
MeV/c, would come entirely from decay of 5 Li, more than 80% of the pairs 

(PL S 1 ^ (PLS) 

would come in the gate v p > Va and this correlation function would be 
much higher than the complementary gate. Since this behavior is not observed, 
it should be safe to conclude that the protons and a-particles contributing to 
the correlations in Fig. [3] (in the region q ~ 54 MeV/c) must have another origin 
than decay of 5 Li. In Ref. PH]> it is assumed that these particles are emitted 
independently from the PLS source, and that the correlations originate from 
final state interaction between particles emitted close in space and time. 
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Figure 4: The experimental velocity-gated pt correlation of Ref. [23! ■ The 
arrows in the figure indicates the positions of different fragment decays (see 
Ref. mi). 



4.3.2 Application to intermediate velocity source emission 

In this section we present numerical simulations similar to those presented in 
section f4. 3. II but now for a different reaction and experimental filter. In Ref. 
PU] velocity-gated small angle two-particle correlation functions were used to 
deduce the emission chronology of protons, deuterons, and tritons, from the 
E/A = 61 MeV 36 Ar + 112 ' 124 Sn collisions. Figure H] illustrates velocity-gated 
p ^correlation functions for the 36 Ar + 124 Sn reaction. In this case the angular 
coverage is 30° < 6 < 114°, and the source of particle emission is assumed to be 
the intermediate velocity source (IS). In particular, we have used the parameters 
summarized in Tabled 

The results of the simulations are presented in Table |31 We see that also 
in this case the velocity condition from the 4 He frame holds, even though for a 



Fragment 


T 


< v > /c 


er„/c 


source 


(MeV) 






TLS 


7 


0.01 


0.01 


IS 


20 


0.18 


0.02 


PLS 


7 


0.31 


0.01 



Table 2: Parameters used in the simulation of section 14.3.21 
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Fragment 
source 


All 
pairs 

(PLS) ^ (PLS) 


30 < 
Ek 

Detected 


Olab < 114° 
> 10 MeV 

(PLS) ^ (PLS) 
Vp > Wa 


TLS 
IS 
PLS 


52.5% 
56.5% 
53.2% 


0.16% 
14.97% 
0.52% 


56.4% 
56.9% 
58.9% 



Table 3: The fraction of pt pairs emitted in the two-body decay of 4 He which 
have v p PLS ^ > v[ PLS \ The last column gives the fraction of "detected" pairs. 



substantial part of the pairs (~ 43%) the triton gets a higher velocity than the 
proton in the IS frame. This means that if the protons and tritons in the region 
q ~ 24 MeV/c, would come entirely from decay of 4 He, then more than 56% 

(IS) (IS) 

of the pairs would come in the gate v p > v{ and this correlation function 
would be higher than the complementary gate. Since this is not observed it 
should also in this case be safe to conclude that the protons and tritons con- 
tributing to the correlations of Fig. 0] (in the region q ~ 24 MeV/c) must have 
a different origin than decay of 4 He. In Ref. [23j it is assumed that these par- 
ticles are emitted independently from the IS source, and that the correlations 
originate from final state interaction between particles emitted close in space 
and time. 



5 Summary 

Two-particle intensity interferometry is an important tool to access information 
on the space-time characteristics of particle emitting sources at intermediate en- 
ergy heavy ion collisions. In particular, when pairs of non-identical particles are 
detected in coincidence, particle- velocity-gated correlation functions can be used 
to establish the emission time sequence, in a model-independent way. Since the 
particle- velocities have to be calculated in the frame of the emitting source, the 
main uncertainty in the method comes from the uncertainty in the source veloc- 
ity, mostly originating from impact parameter averaging implicitly performed 
in experiments. We have demonstrated how the error due to uncertainties in 
the source velocity can be estimated by numerical simulations. A presented ex- 
ample of such an investigation, makes it plausible that for most cases, the error 
will be small, of the order of a few percent. 

A more serious problem may be posed by the fact that for the particle- 
vclocity-gated correlation function method to work, the particles have to be 
emitted independently by the source. Clearly, this is not always the case, as 
particle pairs may originate from two-body decay of excited fragments. We 
have demonstrated that the kinematical signature of two-body decay is strong 
in the particle- velocity-gated correlation functions. This leads to an expected 
behavior, namely that the correlations in the gate where the lighter particle 
has the higher velocity are stronger than in the complementary gate. This 
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ensures that when the contrary behavior is observed, the dominating origin of 
particles must be other than two-body decay. As an example, we have shown two 
experimental cases where it has been possible to deduce the emission chronology, 
even if there is a contribution from two-body decay. 
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